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ABSTRACT 


A model is developed to simulate the public transit system on a corridor 
with Object Oriented Approach This model has the capability for dynamic 
scheduling of vehicles on a major corridor of a metropolitan city This 
model simulates the public transit system on any route network, and 
suggests the optimal utilization of vehicles which minimizes the passenger 
waiting times at all stops This model has a flexibility with respect to arrival 
of passengers to a stop A stop may have any one of Exponential, Normal, 
or Gamma probability distributions of passenger arrival, and different 
parameters of distribution This model displays the location of buses on the 
corridor of a city network at every five minute interval using concurrent 
mode of animation This model supports the interactive graphics which 
enables the user to obtain the information of stops and buses, or to quit the 
simulation process at any instant of time. This model is tested with a case 
study of DIehi ring road, on which the public transit system is simulated A 
tentative timetable covering the morning peak period is prepared with 24 
buses and is fed to the model The model has suggested 19 additional 
buses, and revised the timetable. 
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CHAPTER I 
mXRODUCTION 


1.1 PROLOGUE 

Public tiansport has become an integral and essential feature 
of cities in both developed and developing countiies In populous and developing 
countries like India, where urbanization trend is on the increase, mass transit system 
like commuter rails and bus transportation are mostly unavoidable. In most of the 
Indian cities, public transportation is inadequate and congestion is seveie Despite 
massive infusion of investment, the transportation crisis in metropolitan areas has 
been growing Traffic congestion, air pollution, urban sprawl, and the need for 
energy conservation provide convincing evidence that drastic actions must be taken 
to improve the mass transit system in India in order to alleviate these pioblems 
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In India, among various means of public tianspoitation, the 
bus tianspoit is most impoitant and dominating because of its dooi-to-door 
accessibility and flexibility in opeiation The success of public bus Iranspoit 
howevei depends on the efficient opeiational management looking to the need of 
the passengeis The demand foi bus tianspoitation is veiy high, tup length is shoit, 
and it IS a cheapei mode because of the subsidy piovided by the government in 
diffeient ways. Bus tiansportation system in metiopohtan cities involves the 
movement of laige number of people between lelatively small numbei of given 
locations Bus transpoitation system should piovide adequate mobility to vaiious 
locations to satisfy the demand at all the locations 

In a metiopohtan aiea the bus loutes aie basically designed 
in two diffeient ways direction oriented - wheie the bus moves along majoi 
coriidois in ceitain diiectioiis; destination oriented - wheie the loutes aic designed 
to satisfy the demand fiom one node to the othei In diiection oiiented, the i6?ea is 
that the peisons board the bus at diffeient nodes on majoi coiiidois, and move with 
high speeds in specific diieclions. These loutes generally have heavy demand and 
high fiequency of opeiation The destination oiiented loiiles aie designed to satisfy 
the demand among major 0-D pans and piovide accessibility to various nodes. The 
demand at these major nodes will be geneially high, howevei the demand at 
various intermediate nodes may not be high. Among vaiious issues of diiect 
concern such as loute design, bus allocation, bus scheduling, ciew scheduling, and 
maintenance scheduling^ the design of bus loutes and bus scheduling may be 
consideied as of prime impoitance. 



y 

111 the pieseiit ]jlaiinmg methods, aftei the loiite patli is 
designed, the scheduling plans aie piepaied based on level of seivice, noimally 
defined in teims of load factoi As the daily demand inatiix foi the publi'c tiansport 
IS available, the niimbei of bus tups on a paiticulai loute aie deteimined foi the 
whole day. The numbei of buses that opeiate on paiticulai loute in both directions 
IS decided subject to the opeiational constraints In this way the timetable is 
piepared foi the day wheie we have mmimuni headway during moining and 
evening peak periods, and on othei times the headway is geneially twice the peak 
headway 

1.2 PROBLEM STATEMENT 

Economy in opeiation, leliable, and adequate level of 
service to the users should be the piimaiy objectives of any bus transport 
management It is imperative that the desired level of services to tiie passengers 
should be met by the minimum amount of lesouices. In order to achieve the above 
objectives it is necessary to investigate the following 

• Stiuctuimg of the routes in order to meet the demand in an efficient maniiei on 
existing transit network. 

• Deleimmation of the minimum fleet size for a specified level of service for the 
system. 

• Determination of the schedules with the minimum fleet siZe deteimined above 
to minimize the overall system cost 

Foi the routes along the major corridors, the specific 


problems which come up are - 
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• The significant vaiiation of passenger demand fioni houi to hour 

• The desiied fiequency diiiing the peak period is quite high 

• The arrival of passengeis at the bus stops follows diffeient piobabihty 
distiibutions which is complex to model 

Foi such loutes - the load factoi becomes veiy heavy diuiiig 
ceitain peiiods, passengei waiting time at stops also inci eases consideiably causing 
discomfoit, inconvenience, etc The demand on these loutes is also significantly 
affected by some othei factors such as closing of cinemas, schools, factoiies, 
offices, etc 

In oui piesenl bus system opeiation, a fixed numbei of buses 
aie opeiated as pei a pie-determined time table and no consideration is given to the 
contingencies or bus bieakdowns. In tiaiispoitation and traffic engineering lot of 
electiomc lustiumentation is deployed, but a veiy little information is lepoited to 
improve the public bus tianspoitation system 

To impiove the system of opeiation, it is desiiable to have 
dynamic scheduling based on actual demand pattern. 

1.3 OBJECTIVES OF THE STUDY 

The main objectives of the study includes the following 
tasks 1 elated to the scheduling of public bus tiansportation system on a route 
network of a city. 

1. Geneiation of passengeis at different stops as per diffeient probability 


distiibutions. 
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2 Simulation of the movement of buses fiom one stop to the othei, and scanning 
of passengers in the bus and at the stop 

3 Estimation of statistics foi passengeis with lespect to waiting time 

4 Collection of infoimation about the queues at diffeient stops and decide the 
optimum scheduling of additional buses that may be available 

5 Giaphic display of movement of buses on loiite uetwoik at a fixed interval of 
time, and infomiation of stops and buses 

1.4 LITERATURE REVIEW 

Anderson et al , (1979) have described a flexible simulation 
model of an uiban bus loute m peak hotii tiaffic The simulation piogiain has been 
piesented with paiticulai emphasis on the facility foi mteiactive contiol of the 
simulation run and foi presentation of oveiall and detailed statistics fiom the 
simulation An independent statistical piogiam system based on the mathematical 
models underlying the simulation model has been used to piovide input paianieteis 
for the simulation The piime aioplication of theii simulation model is to lun the 
progiam together with an on-line system that collects data fiom the real bus traffic 
and piesents its status on a display If a difficult situation occuis, the simulation 
piogiam can be consulted foi testing diffeient ways to tackle the pioblem They 
have built the model to handle one bus loute that may consists of seveial service 
valiants 

Fiiedman(1976) foimulated a mathematical model of a 
general public transpoitation netwoik, and suggested an optimized pioceduie to 
schedule the buses’ departuie times. The authoi has assumed that the netwoik 
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opeiates uiidei known deteiministic conditions the lines’ loiites, the stations’ 
locations, the niimbei ofpassengeis and the tiaffic intensity aie given Foi the sake 
of lucidity the model, the model first developed undei the assumption that the 
loading as well as unloading times of passengeis aie zeio Afteiwaids the 
consequences of the relaxation of this supposition aie investigated The 
perfoimance index of the system is the aveiage waiting time of passengeis, and it 
has to be minimized by a piopei choice of decision variables - i e. the buses’ 
depailuie times - without violating the constiaiiils of piesciibed number of buses, 
and diiveis 

Fiiedinan (1975) suggested a lapid solution to the apparently 
simple problem of finding the minimum numbei of buses that aie required to keep 
up a two way lane with piesciibed depaiture and tiavel times with the following 
assumptions 

• The buses tiaverse only the given lane. 

• A bus that ai i ives fust at a tei minal will be the fii st to leave it. 

No empty luns buses fiom one stop to the othei This procedure foi finding the 
minimum numbci of buses that aie needed to maintain a two way lane was utilized 
as pait of a dynamic programming algoiitlim which mean to establish the optimal 
values of the depaituie times of buses without violating the constraint of having a 
given numbei of them. 

Dhingra (1983) has developed a simulation model foi the 
scheduling problem of public bus transportation system, which simulates the flow 
of passengeis on all the nodes of the netwoik and the movement of vehicles on all 
the loiites of the netwoik, theie by taking into account the interaction effects of 
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ovei lapping, ciossing, meiging, and diveiging routes Routes have been designed 
based on liavel demand, and assumed to be independent of schedules. The use of a 
vehicle over multiple loutes has not been coveied. 

Cedei (1986) piovided alternative methods foi constiucting 
bus timetables using passengei load data, and attempted to fulfill six nia|oi 
ob)eclives to evaluate alternative timetables in teinis of lequiied lesoinces, to 
impiove the coiiespondence of bus depaituie times with passengei demand, to 
piovide alternative timetables foi the scheduleis’ use in specific scheduling 
situations, to peimit diiect bus frequency changes foi possible exceptions (known 
to the scheduleis) which do not lely on passenger demand data, to allow the 
constiuction of the timetables with headway smoothing techniques; and to integrate 
diffeient headway setting, and diffeient timetable constiuction methods. 

Dumas et ah, (1990) piesented an algoiithm that solves the 
pioblem of finding the vehicle schedule which minimizes total inconveniences foi 
tiavcl along a fixed path, wheie seivice tames at nodes aie constiained by time 
windows, and wheie inconvenience is modeled using convex functions of the 
seivice times. It has been shown that the complexity of the algoiithm, expiessed as 
a numbei of unidimensional minimizations, is on the oicler of the numbei of nodes 
foi convex inconvenience functions. An extension was pioposed which includes 
lineal costs on waiting tunes, and disciete service time vaiiables 

Scheele (1980) foiimilated a pioblem to determine the 
complete travel pattern and to decide which bus frequencies to use on the various 
lines The pioblem is formulated as a non-linear progiamming pioblem in the foim 
of a compound minimization pioblem. The assumptions made to simplify this 



pioblem aie - the netwoik of stieels on which ceitam bus lines have been set up 
IS given, the total numbei of buses is given, the total demand foi bus ti anspoi tation 
IS given in the foim of the marginal totals of an oiigin-destination maliix An 
ileiative algoiithm to solve this problem is developed 



CHAPTER II 


SIMULATION MODEL 


I 

2.1 USE OF DIGITAL SIMULATION 

The mass transportation network of a major urban aiea 
necessarily involve a large number of logical sequence of decisions in the flow logic 
The scheduling process being intricate requires the use of digital simulation Digital 
simulation may be defined as, ‘a numerical technique for conducting experiments on 
a digital computei with certain types of mathematical and logical models that 
describe the behavior of the complex system ( oj’ some part thereof) over 
periods of real time’ ( Naylor et al, 1967 ) Simulation serves different functions in 
system analysis and design Digital simulation has the following additional 
advantages - 

1 Data can be generated which enables the researcher to study and experiment 
with complex internal interactions of a given system and thus bettei 
understanding of the nature of interactions. 
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2 Simulation of complex systems yield valuable information for identifying the 
dominant parameters governing the system operation 

3 It IS useful in validating the component models in conjunction with the global 
modeling of a system 

2.2 OBJECT ORIENTED PROGRAMMING 
(OOP) APPROACH. 

This model is developed with OOP approach in C on 
Turbo C platform This section gives a brief idea, of OOP appioach, and the use 
of this approach m this model This section helps in undei standing the OOP 
approach, and thus the piogram veiy easily, however, the readers having the 
sufficient knowledge of OOP approach may skip this section with out any 
difficulties 

OOP IS a way of organizing programs The emphasis is on the way 
progiams are designed, not on the details of individual operators, In particular, 
OOPs are organized around objects which contain both data, and functions the act 
on that data A class is a template for a number of objects The fundamental idea 
behind OOP approach is to combine both data, and the functions that operate on 
that data into a single unit, which is called an object An object’s functions called 
member functions in C , provide the only way to access its data If for example, 
you want to read a data item in an object, you call a membei function in the.object 
It will read the item, and return the value to you The data is hidden, so it is safe 
from accidental alterations. Data, and its functions are said to be encapsulated into a 
single entity If you want to modify the data in an object, you know exactly what 



functions interact with it the member functions in the object No other function can 
access the data This simplifies writing, debugging, and maintaining the program 

I 

A C progiam typically consists of a number of objects 
which communicate with each other by calling one others member functions 
Thinking in terms of objects, rather than functions has a surprisingly helpful effect 
on how easily programs can be designed This results fiom the close match between 
objects in the programming sense, and objects m the real world ( Lafore, 1994 ) 

2.3 MODEL DESCRIPTION 

In the previous section the concept of OOP is briefly 
introduced This section, now describes the model, and discuss how the concept of 
OOP is used to develop this model 

This model covers the simulation of public Bus 
transportation system on a route network of a city given travel demand using 
heuristic approach The stochastic nature of flow of passengeis and vehicles in a 
network necessitates the use of simulation modeling to achieve the objectives 

# * I 

mentioned in section 1 3 The simulation model developed in this study simulates at 
micro level the ariival of passengers at all Stops as per different probability 
distributions, and the movement of Buses from one Stop to the other, The salient 
featuies of the model that accounts various complexities of the system are as 


follows - 



15 


1 The model simulates simultaneously the flow of passengers at all the Stops of 
the network, and the movement of flow of vehicles between all Stops of the 
network 

2 The model takes into account the stochastic nature of passengei ariivals For 
each node, a probability distribution and its parameters aie to be specified 

t 

reflecting the ai rival gap characteristics Exponential, Normal, and Gamma 
distributions have been associated with each Stop for the complete simulation 
run 

3 This model suggests the schedule of an extra Btis from any of the terminals 
which may have exti&Bitses, based on some heuiistic approaches Also the user 
can accept or reject the suggestion made by this simulation model 

4 This model, for every five minute interval displays graphically the locations of 

« > 

Buses on the network, and the information of Stops and Buses using concurrent 
mode of animation 

5 The graphic display is highly usei interactive, that the usei can Stop the 
simulation run at any instant of time, and quit the graphic display 

The progiam developed for this model is divided into a 
numbei of objects Each Bus Stop on the network is considered as one object of 
class Stop. Similaily, each Bus is also treated as one object of class Bus This 
program can handle a system with any number of Stops and Buses. If, foi example 
there aie M number of Stops, and N number of Buses, this progiam dynamically 
creates memory for M Stop objects, and N Bus objects The data structures 
associated with Stop objects. Bus objects are Linked Lists, and Queues, to maintain 
the passengers in each Stop, and each Bus, and to maintain the queue of departing 



Buses at each Stop The detailed cfescnption of Stop objects, and Bus objects will 
be given in later seetions 


2.4 MODEL INPUT AND OUTPUT 

This model takes input fiom four different files - 

1 d 

1 In first input file will have the 'information of all Stops on the route network 
The fust element in input file, 'STOP_XY.IN' will be MAX STOP , the 
maximum numbei of Stops under consideiation on the network of a city The 
information of each Stops with respect to its sequence number, x and y 
coordinates, the destied probability distribution of passengers arriving at that 
Stop ( EX for Exponential, NR foi Normal, and GM foi Gamma) and its 
paiameters, and the total number of Buses available at that Stop have to be 
stored m tins file 

2 In the second input file the information of all Buses will be stored The first 
element in the input file ‘BUSJfNFO.IN’ will be B(/S_CAP, the maximum 
capacity of Bus The information of the each Bus (that are presently being 
operated in the system) with respect to its sequence number, staring terminal, 
starting time, and direction of movement on the network ( CW foi clock wise, 
and ACfV for anti clock wise ) will have to be stored in this file, 

3 , The input file No 3, ’TIMEJTAB.IN’ will have the departure times on all Buses 
at each of the Stop on the network of the city The Bus id number, and its 
departure time, following the id number of Stop will be the right format In this 
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way the departure times of all Buses at each Slop will have to be stored in this 
file 

4 The input file No 4, ‘IJNKTIME IN ’consists of travel times on all the links of 
the netwoik The link in terms of starting Slop, and ending Slop has to be 
mentioned in the input file, then follows the tiavel time on that link, and the 
dll ection of that link, either CW or ACW 

5 The input file No 5, ‘OD_MATIN ’ will hold the OD matrix of size 
MAX STOP in which each element represents daily demand of passengers from 
one Slop to the other 

The output of this simulation model is in the foiin of giaphic 
display mostly, and can be classified into thiee categoiies - 

t 

1 At every five minute interval the graphic display of the network of the city, and 
the locations of Buses on the network 

2 At every five minute interval the graphic display of Stop infoimation with 

lespect to the sequence icf number of the Slop, the number of passengeis 

waiting at that Stop, the average waiting time of those passengers who board 

the last departed Bus from that Slop, the departure time and sequence numbei 

of the next Bus fiom that Stop, and the number of spaie Buses that are available 

1 •* 

at that Slop 

3 At eveiy five minute interval the giaphic display of Bus information with respect 
to the sequence id number of the Bus, the no of passengers occupied the Bus, 
the departure time and sequence number of previous Stop, and next Stop 

4 Graphic display of suggestion of additional Bus that can be scheduled from any 
of the teiminals if satisfies all the conditions 



' ll 

5 The final output file gives the schedule of all Buses including the newly 
scheduled additional Bt/ses, from all Stops 

2.5 MODEL FORMULATION 

This section explains in detail the formulation of the model, 
and Its development As explained in section 2 2, depending on the size of the 
problem, this model creates the Stop objects, and Bus objects This section is 
divided into seven parts for the bettei understanding 

Tn part one the general flow diagram is illustrated which 
gives the idea of how the diffeieiit modules or functions of the model are intei- 
related In paits two and three of this section the data and member functions 
associated with each Stop object, and Bus object are explained In pait four the 
beginning of simulation process is described In part five it is explained how the 
passengers will be generated at each Stop in the model. In part six the techniques of 
scanning the Stops, and Buses is discussed In part seven the stiategy to fihd the 
best Stop from which an additional Btts could be scheduled is descnbed In part 
eight the method used in this model to compute the location of Buses on the route 
network is explained The last part i e part nine of this section discusses the 
teimination critenon of the model 

2.5.1 GENERAL FLOW DIAGRAM OF MODEL 

The rough flow diagram shown below gives a precise idea how the 
model works, or how the various fimctions are inter-related to each other However 
each function and its sub-functions will be explained in latei sections with the help 


of flow charts 
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Fig. 2.1 General Flow Diagram of the Model 
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2.5,2 STOP OBJECTS 

There will beMAXSTOF number of Stop objects created by 
the piogiam Each Siop object will have the following data items, and member 
fliiictions to access the data from outside functions, and to inteiact with its data 
a) DATA ITEMS : 

• id Sequence number of the Stop 

• x co, y CO X, and Y cooidinates of Stop 

• dfstf Desired probability distribution of passenger arrivals at Stop 

Can be any one of Zfxponential, jVoimal, oi Gamma 

• part First parameter of the desired probability distribution 

• par2 , Second parametei of the desired probability distribution 

• tot_bus . Total number of Busts a Stop has 

• * , 

It is equal to'zeio if a Stop is not a terminal 

• free bus . Spare or additional 5w5es that are available fionia5?o/i 

• tot_wait . Total waiting tune of all passengers served at the Stop 

• civg_\mH . Avei age waiting time of those passengers who board the 

last Bus departed from Stop at any instant of time 

• PassQue The queue of passengers arrived at the Stop as per 

paiticulai distribution Each element of PassQue will consists of the arrival 
time of a passenger, and his/her destination When ever a passenger is 
generated, he/she will be stored in ibe PassQue of lespective Stop, and 
similarly when a Bus departs the Stop he/she will be deleted from the 

PassQue, and stored in PassLmk of that Bus. 
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BusQue The queue of departing Buses from the Stop. -in the 

/ 

increasing order of their departure times Each element in BusQ}te represents 
the Bns id numbei, its departure time, and its direction of movement 
In other way, this is the time table for the Stop Whenever a Bus 
departs the Stop^ the first Bus in the BusQue will be deleted, so that the 
next Bus will be ready to depart 
b) MEMBER FUNCTIONS ; 

• gel_data( stop no, x_co, y_co, distr, pari, pai 2, tot_bus ) 

This function stores the data . obtained from the input file No 1 into the 
1 espective Stop object 

• getjinietah( bus no, busjime, dim ) 

This function stores the details of departing Bus te the Bus number, its 
depai tut e time , and its diiection of movement in the BusQue of object Stop 

• retum^idO 

This function returns the id number of a Stop, and is useful to check the id 
numbei of a desired Stop 

• return J'trst_bus( ) 

This function returns the id numbei of Bus which is next to depart from the 
Stop. 

• letiirn J'irstJime() ■ 

This function returns the departure time of Bus which is next to depart from the 
Stop. 

• return Jasi_time( ) , 



This function i etui ns the departure time of last Bus from the Stop This function 
IS useful to check whether the Stop is empty or not 
move Ji! st_bus( ) . 

This fijnction deletes the first Bus from BvsQue when it departs the Stop, i e 

f 

when ever a Bus departs a Stop, BusQue of that Stop will be updated 
add j)ass(} ' 

This function stores the randomly geneiated passenger to the PassQue of the 
Stop 

delj)ass( ) 

This flinction deletes the passenger from the PassQue of the Stop when the 
passenger boards a departing Bus 
populatiou( ) 

This function returns the population of passengers, i e the number of 
passenger s in tiie PassQice of the Stop, 
return Jot jvait( ) 

this function returns the toijvait of the Stop 
retuinjivgjvtime{ ). 

This function i etui ns the avgjvajt of those passengers who boaid the previous 
departing Bus from the Stop 

§ 

givejreejus( ) 

This function reiwrns free J)us, the number of spare buses that the Stop has 
decrement free J3vs( ) . 

This flinction decrements the free bus of the Stop, when an additional Bus has 
been scheduled from the Stop 



• retii)V_x_co( ) 

This function returns the x_co of the Stop 

• retmv_y_co( ) 

This function returns the y_co of the Stop 

• ielit}njno_btis() 

This function returns the number of Buses that are yet to pass through the Stop 
This fimction, when called checks the BiisQue of a x^lop and returns the number 
of Buses in BusQue. If theie aie no Buses, i e if the last Bus has departed, this 
function 1 etui ns zero which implies the Stop is empty 

• jefuin 

This fi-inction returns the first parametei of the probability distribution at the 
Stop 

• return j>ai 2 () 

this function leturns the second parameter of the probability distribution at the 
Stop 

• check_distr(d ) : 

This function compares the distribution d given to this function as argument 
with the desiied distiibution at the Stop, and leturns YES if both are same, and 
1 etui ns NO if not 

2.5.3 BUS OBJECTS 

MAX_BUS, the numbei of Bvs objects to be created by the 
piogiam, depends on the loi bus of each Stop 'The MAX BUS is nothing ’but the 
sum of iot_biis of all Stops. But FLEET_SIZE, the numbei of Buses being operated 
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on the route network, is initially equal to the number of Bus objects given in the 
input file No 2 FLEET _SIZE will be incremented by one unit when an additional 
Bus IS scheduled The value of FLEETJUZE can not exceed -as no 

f 

additional Bus can be scheduled if there aie no free _bus at each Stop Now let us 
look into the data items, and member functions of each Bus object 
a) DATA ITEMS ; 

• id Sequence number of a Bus 

• tei mjd' Sequence number of starting tei minal of the Bus 

• dun. Direction of the movement on the route network 

rif- Clock Wise 01 ACW- Anti Clock Wise 

• prev_stop Sequence nutnbei of previous Stop 

• prevjime Departure tune at previous Stop 

• PassList: The Linked List of passengers who had occupied the Bus 

Each element of PossList compiises the destination of a passenger When ever 

the Bus 1 eaches a Stop, the PassList in the Bus will be scanned, and all those 

passengers whose destination is that Stop will be deleted from the PassList, and 

the passengers from StopQue will be added to the PassList of the Bus subjected 

to the capacity constraint of the Bus 

0 

• StopQue: The queue of Stops to which the Bus has to travel Each 

element of StopQue comprises the id number of a Stop, and the departure 
time of a Bus at that Stop, in the increasing ordei of time In othei way this is 
the time table for the Bus When ever a Bus reaches a Stop, the first element 


will be removed from the StopQue 



MEMBER FUNCTIONS ; 


get_data( id, term xd, dirn ) 

This function gets the data, the sequence number of a Bus^ its starting terminal, 
and its direction of movement from input file No 2 and stores them into the Bus 
object 

getjmtetable( stop_\6, depjime) 

This function stores the data, the sequence number of Stop, and departure time 
at that Stop into StopQue of the Bus 
returnj^( ), 

This function retui ns the id number of the Bus, and useful to check the /d of a 
Bus 

del_bus_siop( ) 

This flinction deletes the fust Stop from the StopQue of a Bus when the Bus 
departs that Stop, and thus updates the StopQue of that Stop 
1 etnrn_stop_vd( ) 

This function returns id of the first Stop in StopQue of the Bus, i e the 
approaching Stop number at which the Bus is going to Stop 
reUini_stopjime( ). 

This function returns the departure time of the Bus at the coming Stop 
add_pass( destination) 

This function adds one passenger from PassQiie of a Stop to PassLisi- of the 

i 

Bus 


del j)ass(). 
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When the Bus reaches a Stop, this flmction will be called This function deletes 

** , 

all the passenger from PassLtsi pf the Bus whose destination is that Stop 

• population( ) 

This function, when called returns the population of PassList of the Bus, i e the 
number of passengers occupied the Bus 

• stole j)revjiata( ) 

This function stores the previous Stop id number, and the departure time of the 
Bus 

• leturn _prev_stop(). 

This flmction returns the sequence numbei of a previous Stop 

• return _prev_Ume( ) 

This Rinction returns the departure time of the Bus at the previous Stop 

• return _no_Stops( ) 

This function returns the number of Stops the Bus has to pass through. In other 
view, this funcion gives the number of Stops in StopQue of a Bus 

• check_dir(d ) 

This fianction compares the direction of Bus and the direction d passed to the 
flmction as an argument and returns YES if they aie same, and NO if not 

2.5.4 BEGINNING THE SIMULATION 

Once the Stop and Bus objects are created, and the input is 
given to the model, the model calculates the systemjme, the time to start the 
simulation process, In this model, the system time is defined as the time 30 minutes 



befoie the minimum departure time of first Bus of all Stops Expressing 
mathematically, 

system time = Min { dep_time of//r5/_^>?«' from iS'/o/'[i], V i }- 30 minutes 

Hence, at system time, the model starts simulating the public 
transit system on a given route network Thereupon, system June will be 
incremented by 5 minutes During each of this five minute interval, passengers will 
be generated at all Stops,, which will be explained in detail m section 2 5 4 each 
Stop and Bus will be scanned at every five minute interval and boarding and 
alighting of passengers will be earned out The section 2 5 5 explains the same 
process 

Also, each Stop will be examined during each five minute 

interval to see whether an additional Bus can be scheduled from that Stop oimot If 

/ 

theie are more than one Stop fiom which an additional Bus could be scheduled, then 
that Stop which has the maximum saving in waiting time will be chosen to provide 
the additional Bus The evaluation of alternate strategies to find best stop fi'om 
which an additional Bus could be scheduled will be explained in section 2 5 6 At 
every increment of system time, the locations of all Busqs, on the network will be 
displayed graphically To display the Busq& on the network the position of each Bus 
at systemjtime is to be calculated The section 2 5 7 briefly explains the same 
Accordingly the model will carry through the simulation process until the 
termination criterion is found, which will be explained in section 2,5,8 
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2.5.5 GENERATION OF PASSENGERS 

This section explains how the passengers are generated 
randomly at each Stop as per specific probability distributions and parameteis 
obtained from input file No I This model can handle the passenger ariivals as pei 
Exponential, Normal, and Gamma distributions Each Stop can have different, or 
same probability distributions, and the parameters In input file No 1 , the desired 
pi obability distributions, and its parameters have to be specified foi each Stop 

The function generate jpass( ) generates the passengers, 
assigns them a random destination, and stores them in PassQue of respective Stop 
At each five minute increment to systentjww, this function will be called for each 
Stop if there is at least one departing Bm from that Stop within 30 minutes from 
system Jime In other woids, passengers will not be generated at any Stop if there is 
no Bus to depart that Stop within half an houi fiom system Jtme This model 
assumes that passengers start ariiving as per pre-decided probability distributions at 
any Stop thirty minutes pnoi to the departure time of Bus from that Stop This 

^ ^ i 

implies that all passengeis are aware of Bus timings fiom each Stop Now let us 
have some /t/ea on how the passengeis are genet ated in this model as per different 
probability distributions at Siop& Though the flow charts will be used in later 
section to understand the logic in a better way 
a ) EXPONENTIAL DISTRIBUTION : 

The function generate _pass( ) generates the passengers at 
any Stop as per Exponential distribution if ‘EX ’ is entered in input file No 1 as 

** I 

distribution at that Stop Exponential distribution will have mean as its parameter. 
The parameter X, the mean or inter-arrival rate should be greater than 0, in any 



case If the value of A entered is less than zero, this function takes a 2 0 as the 
default value of A , and proceeds with the flirthei calculations In this function 
inverse-transfoim technique is used to geneiate Exponential random variates 

Let X have the exponential distiibution with mean X The 
distribution function is - 

F(x) = ifx>0, 

= 0 otherwise 

So to find F"\ we set u = F(x), and solve for x to obtain 
F-'(u) = - X In(l-u) 

Thus to geneiate the desiied random variate we first geneiate a U ~ U(0, 1), and 

“•* I 

then let x = -X In U It is possible in this case to use U instead of 1 - U, since 1 - 
U and U have the same U(0,1) distiibution. This saves a subtraction (Law, and 
Kelton, 1991), 

b ) NORMAL DISTRIBUTION : 

The fiinction generate_pass( ) generates the passengers as 
per Normal distribution if NR is the desired distribution at any Stop The two 
parameters of Noimal distribution, i e mean p, and variance have to be entered in 
the input file This function howevei assumes the minimum values of mean and 
variance as 0 and 1 respectively 

The polar method is used in this function to generate N(0,1) 
random variates (Law and Kelton, 1991) Which is as follows - 
1 Generate U i , and U 2 as ITD U(0, 1 ) Let V, = 2U, - 1 for i = 1 ,2 


LetW- Vi^ + V 2 ^. 
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2. If W > 1, go back to step 1. 

Otherwise, let Y = V[(-21nW)A¥]. Xi=ViY, i=l,2. 

Then Xi, i = 1,2 are lED N(0,1) random variates. 

Having X ~ N(0, 1) , we can obtain X' ~ N (n,a^) by setting X' = p + a X. 
c ) GAMMA DISTRIBUTION : 

The function generate_pass( ) generates passengers as per 
Gamma random variates if the desired distribution at a Stop is entered as GM. The 
parameters of Gamma distribution are: shape parameter a >0, the scale parameter P 
>0. The mean and variance of Gamma distribution are aP, and aP^ respectively. 
This functions takes 1, and 1 as the minimum values of a, and p. The algorithm, 
recommended by Cheng is used in this function to generate Gamma random 
variates, which generates gamma(a,l) random variates ( Law, and Kelton, 1991). 

The prespecified constants are a = 1 / V(2a - 1), b = a - In 
4, q = a + 1/a ; 0 = 4.5 ; and d = 1 + In 0. 

1 . Generate Ui, and U 2 as IID U(0, 1). 

2. Let V = a ln[U,/ (I-U 2 )], Y = ae'', Z = Ui'Uz . and W = b + qV -Y. 

3. IfW + d- 0Z>O, return X = Y. Otherwise, proceed to step 4. 

4. If W > In Z, return X = Y. Otherwise, go to Step 1 . 

Given X ~ Gamma(a, 1) random variate one can obtain X' ~ 
Gamma(a,p) variate for any P > 0, from X' = PX. 

ASSIGNING THE DESTINATIONS : 

0 

As stated earlier, the function generatejpass( ) generates the 
passengers as per pre-decided distributions, and assigns them a destination 



34 


randomly. The destinations will be assigned to passengers with the help of OD 
matrix entered in input file No. 5. Each element of OD matrix represents the daily 
demand from one Stop to the other. A function OD_manipulation( ) is used to 
manipulate the initial OD matrix so that each element of manipulated OD matrix 
represents the cumulative percent of trips (daily demand) from one Stop to the 
other. 

In the function generate_pass( ) the manipulated OD matrix 
will be used. A random variate U ~ U(0,1) will be generated and if that random 
number matches with any block of cumulative percent of trips of all destinations for 
a given Stop, that destination will be the passenger’s destination. The flow chart will 
illustrated in section 3.3.4 to explain the same logic. 

2.5.6 SCANNING STOPS AND BUSES 

In this section let us see how this model scans each Stop, and 
each Bus, and how it alights/boards the passengers ffom/to each Bus. This process 
will be repeated for each Stop, and for each increment of system Jime. 

th 

For each increment of system Jime, following steps will be taken For i Stop, 

1 . Calling return ) of Stop\i\, the id. of that Stop is stored in Stopjio. , 

r 

2. The passengers waiting at Stop[i\ is obtained by calling Stop[\\.population { ) 

3. From the BusQue of the departure time of Bus , and its id. are noted as 

bus Jime, and busno respectively. This implies, the Bus whose id. is busjio 
departs the Stop[i\ at bus Jime. Now the system Jime is compared with 


bus time. 
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4. The system time equal to busjime implies that that Bus is about to depart. 
Now the boarding and alighting of passengers to and from the Bus at Stop\i\ 
will be done. 

5. The system time greater than bus time implies that that Bus has already 
departed the Stop. In this case also, the alighting and boarding of passengers 
from and to the Bus at Stop[\\ will be carried out. 

6. The system Jime less than the busjime implies that the Bus has not arrived the 
Stop]}]. In this case, nothing will be done as there is no Bus at Stop\\\. Now the 
same process from step-1 will be repeated for next Stop, {.&.( i +l/^ Stop. 

BOARDING AND ALIGHTING OF PASSENGERS : 

Alighting and boarding of passengers from and to a Bus will 

be done only when the system time is greater than or equal to busjime. The 

following steps will be taken to achieve this ; 

0 

1 . Scan the passengers in the PassLink of the Bus. 

If there are any passengers whose destination is equal to Stop jio, unload them 
from PassLink of the Bus by calling the member function del j)ass{ ) of that 
Bus. 

2. Call the member function population{ ) of the Bus, obtain the population in that 
Bus, i.e. the number of passengers occupied the Bus. 

3. Compare the population in the Bus with BUSjOAP, the maximum capacity of 
the Bus. If the population is equal to BUSJSAP, then the remaining steps will 
not be executed, as the Bus can not accommodate any ,more passengers. In this 
case, the allow the Bus to depart the Stop. 
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4. If the population in the Bus is less than BUS CAP, load the Bus with the 

f 

passengers from Stop\i\. Check the arrival time of first passenger in the 
PassOue. 

5. Scan the population at the Stop, i.e. the passengers waiting at the Stop\\\ by 
calling the member function population^ ) of Stop\i\. 

6. If the population at Stop\\\ is greater than zero, and if the passenger’s arrival 
time is before the departure time of the Bus, load the passenger into the Bus by 
calling the member function add jpass( ) of that Bus. If the arrival time 
passenger is greater than busjime, allow the Bus to depart, as the passengers in 
the Stop arrived after the departure time of the Bus. After loading the passenger 
to the Bus, delete that passenger fi-om PassOue of Stop\i\ by calling del jpass { ) 
of Stop\i\. 

7. If the population at Stoplyl is equal to zero, i.e. if there are no passengers at the 
Stopli], allow the Bus to depart, as there are no passengers waiting at Stop[\\. 

MOVING THE BUS FROM STOP : 

If there are no passengers at the Stop{\\, or if there is no 
space in the Bus for passengers, of if the alighting and boarding of passengers is 
finished, then allow the Bus to depart the Stop. Call move Jirst_bus{ ) of Stop{\\, 
which deletes that Bus from the BusQue of Stop\\\, and updates the BusQue. 

Once the Bus is moved out of Stop{f\, then the same 
procedure will be repeated for next Stop, i.e. Stop{\ + 1]. If this procedure is 
repeated for all Stop%, then the next step to be under taken is fihiding the best_stop 
from which an additional Bus could be scheduled. This method is described in the 


following section. 
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2.5.7 FINDING THE BEST STRATEGY TO PROVIDE 
AN ADDITIONAL BUS 

This section explains the evaluation criterion to find the best 
strategy to provide an extra Bus in to the route network. This evaluation will be 
carried out at each five minute interval of systemjime. The function 
find_best stop( ) will be called each time, which examines all the Stops to find the 
best stop which can provide the maximum saving in waiting time for the 
passengers. 

SELECTING THE ALTERNATIVES : 

First the alternative Stops for best stop will be traced out based on some heuristic 
conditions, then the best_stop, the best of all alternatives will be found, .and an 
additional Bus will be suggested by the model. The three conditions that have to be 
satisfied for a Stop to be the one of alternative are - 

• The basic condition is that the Stop must have at least one extra Bus. A Stop 
with out an extra Bus can not be considered as an alternative as it is not possible 
to schedule any Bus from that Stop. 

• The condition favoring system operators is that there must be at least some 

minimum number of passengers waiting at a Stop which satisfies the first 

** 1 

condition. This condition is imposed keeping in view the operators’ benefits. If 
there are a few number of people waiting at the Stop will not make any good for 
the operators. Any number which suits the system can be considered as the 
minimum number of passengers. However, in this model, 20 is considered as the 
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minimum number. More clearly, only that Stop which has minimum 20 number 
of passengers waiting at system time is said to be satisfy the second condition. 

• The third condition is that there must be no Bus at a Stop with departure time 
lesser than or equal to system time plus some minimum time period. If a Stop 
satisfies the first two conditions, and only if there is no Bus from that Stop to 
depart within some minimum period of time from system Jime, that Stop can be 
considered as the alternative. This minimum time period can be any thing from 0 
minute, however minimum time'period considered in this model is 5 minute. So, 
a Stop to be one of the alternatives have to satisfy the first two conditions, and 
there must be no Bus leaving that Stop with departure time falling in the five 
minute period from systemjime onwards. 

If there are no Stop% which satisfy all the three conditions, 

then there would be no alternatives, and no additional Bus will be suggested. If 

there is only one Stop satisfying all the three conditions, then no evaluation of 

alternatives is necessary, as there is only one alternative, and that Stop .will be 

$ 

considered as the best_stop, and an additional Bus will be suggested from that Stop. 
Only if there are more than one Stop satisfying all the three conditions, then 
evaluation of alternatives will be made to find the best_stop as explained here. 
EVALUATING THE BEST ALTERNATIVE : 

Evaluation of alternatives will be made only if there are more 
than one Stop which is considered as an alternative, i.e. which satisfy all three 
conditions stated above. For each alternative, first the saving in waiting time will be 
calculated. Then the alternative which has the maximum saving in waiting time will 
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be considered as the best stop. The following steps will be followed to calculate 
the saving in waiting time for each alternative, say for Stop\\\ - 

1. Get bus time, the departure time of next Bus from Stop]\\ by -palling 

t 

return Jirst_time { ) of Stop[i]. Calculate the timejiiff, the difference between 
the system Jime, and busjime. 

2. timejiiff = system Jime - busjime. 

This implies that if an additional Bus is scheduled from Stop\\\, it will save the 
waiting time of each passenger by timejiiffmmxAQs. 

3 . Now get the population of Stop{\\ from population ^ ) of that Stop. 

4. Saving in waiting time at Stop{\\ is the product of population at Stop\jL\ and 
timejiiff. 

5. Similarly get the saving in waiting time at all the other Stops, if the Bus from 

Stop[i] is scheduled at system Jime. Saving in waiting time at each of the Stops, 

other than Stop\i] is the computed as the product of population of that Stop, 
and timejiiff. 

6. Now calculate the total saving in waiting time which is the sum of saving in 
waiting time of all Stops. 

7. In case of movement of Buses on both directions, calculate the total saying in 
waiting time at a Stop in both directions separately, and the direction in which 
total saving is maximum should be taken into account. 

Once all the Stops are evaluated with respect to the total 
saving in waiting time, the bestjtop will be that Stop with maximum total saving in 
waiting time. This model evaluate the bestjtop considering the both directions of 
movement. Thus and so the function ^nd_bestjtop{ ) finds the bestjtop, and the 
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direction of Bus to be scheduled from that Stop. The model displays a message 
graphically which states that an additional Bus is scheduled from the bestjstop at 
the system time., and gives the user the choice to accept it or reject. If the user 
accepts the suggestion, then the new Bus will join the fleet on the network. In this 
case the FLEET_SIZE will be incremented by one, and the jree_bus of the 
best stop will be decremented by one, and the simulation continues further. If the 
user rejects the suggestion made by the model, then the new Bus will not be 
scheduled and the simulation continues further. 

2.5.8 COMPUTING BUS LOCATIONS 

This section explains how the at every increment of 

system time the Bus positions are calculated to display their locations on the 

network. The function getjbusJocations{ ) will be called after the function 

find_best_stop{ ). The following steps will be followed by the fiinction 

*'* . 

get_bus_locations( ). 

At any system time., for every Bus - 

1. For a Bus, from StopQue get the next Stop number, and departure time at that 
Stop using functions return _stopj.d{^ ), and return_stopJime { ). 

2. Get the previous Stop details for a Bus from return _prev_stop{ ), and 
return jprev_time(). 

3. The system Jime equal to prevjime -1, or prevjime implies that the Bus is at 
prevjstop at systemjime. Thus the position of that Bus is that of the preyjstop. 

t 

Now plot the position of^M^ on the network. And go to step 1, repeat for next 


Bus. 
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4. The system time equal to nextjime implies that the Bus is at next_stop at 
system time. So the position of the Bus is that of nextjstop. Plot the position of 
the Bus on the network, and go to step 1 and repeat for next Bus. 

5. The system Jime greater than prev_time, and less than nextjime implies that 
the Bus is in between the two Stops. In this case it is necessary to do some 
calculations to get the position of the Bus. Proceed to step 6. 

6. Compute the time lag, the difference of prevjime, and system Jime. Compute 
the tot dist, the distance between nextjtop, and prevjtop in terms of pixels. 
Obtain the travel time from prev stop to nextjtop from input file No. 4. 

7. Now calculate the desired speed of the Bus, which is equal to totjiist divided 
by the travel time. Get movejiist, distance (in terms of pixels) to move the Bus 
from prev stop is equal to the desired speed of the Bus multiplied by the 
time Jag. 

8. Calculate the slope of Bus movement from prevjtop to nextjtop depending on 
the direction of the Bus. Obtain the angle 9, which is inverse tangent of slope. 

9. Now, at system Jime the current position of the Bus on the network is obtained 
from - 

Bus.xjo = prevjtop.xjo + movejiist x cos 0' 

Bus. y JO = prevjtop.y co + movejiist x sin 0' 

where 0' = 0 if the slope is south east, or north west. 

= 90-0 if the slope is south west, or north east 

10. Plot the Bus position on the network. Go to step 1, and repeat for the next Bus, 
until all Busqs are plotted. 
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2 . 5.9 TERMINATION CRITERION 

The model will proceed with the simulation till the 
termination criterion is found. If the last Bus from a Stop has departed, then that 
Stop is said to be empty. So if all Stops are found empty then the termination 
criterion is supposed to be satisfied. To check whether a Stop is empty or not call 
the member function return_no_5M5( ), which returns the number of Busts in 
BusOue of that Stop. If the termination criterion is found, the model displays a 
message graphically that all Stops are empty, and asks the user to quit the model as 
no further simulation is possible. 
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CHAPTER III 

PROGRAM STRUCTURE 


3.1 ABOUT TURBO C++ 

As Stated in section 2.2 this program is 
developed with OOP approach using and run on Turbo compiler Version 2.0 

in an MS-DOS environment. Turbo C provides a powerful environment called 

0 

integrated development environment (IDE) for creating and executing a program. 
The IDE is completely menu-driven and allows the user to create, edit, compile, and 
run programs using what are known as dialogue boxes. Turbo C requires MS-DOS 

2.0 or higher, and runs on PC XTs, ATs, PC/2S, and compatibles. It needs at least 
640 KB of memory, a hard disk, and a floppy disk drive. The hard disk should have 
some what more than 6 MB of free space. Turbo will also run in Microsoft 

Window DOS box. Turbo works with any 80 column monitor either character 

0 

based ( which usually means monochrome ) or graphic-based ( usually color ). Turbo 
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C can also accommodate a mouse, how ever it is not necessary; one can perform 
all operations from the key board. 

3.2 MAIN PROGRAM 

f 

The main program is divided into 6 major functions, a few of 
which are further divided into sub functions. This section explains the complete 
program with the help of flow charts. The six functions which will be explained in 
later sections are - 
/. graphics display ( ) 

2. create stops( ) 

3. create huses( ) 

4. assign _time_table( ) 

5. odjnaniptdation( ) 

6. do operation( ) 

3.2.1 GRAPHICS JDISPLAYi ) 

The main program of this model calls five functions stated 
above. First function to be called by the main program is graphics _displc^( ). This 
function in itself calls initgraph( ), a library function which must be executed before 
any other graphic-mode function used. This library function switches the computer’s 
display system into appropriate graphics mode. This function requires the ghaphics.h 
header file. After initiating graphics, graphics _display( ) function calls other graphic 
functions used in the program to create windows, push-buttons, and text. The steps 
which will be executed in this function are given below. 


1. initgraph() 
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2. display_main_window( ) 

3 . display_time_window( ) 

4. display_push_buttons( ) 

5. display_text( ) 

3.2.2 CREATE_STOPS( ) 

This is the function that will be called after the 
graphics _display( ). This function 'reads the first element MAX_STOP of input file 
No. 1, and creates the MAXJSTOP number of Stop objects, and then stores the data 
from input file No. 1 by calling the member function get_data( ) of each Stop. This 
function calculates MAXJBUS, the maximum number of Buses that can be used in 
the system which is the sum of tot bus of all Stops. 

3.2.3 CREATEJBUSESi ) 

This function opens the input file No. 2 and first reads 
BUSjCAP, the maximum capacity of a Bus. Then it creates the MAXJBUS number 
of Bus objects. Then it reads the input file No. 2, and computes FLEET_SIZE, the 
number of buses being operated initially, which is the number of Buses defined in the 
input file No. 2. Then it calls get_data{ ) of each Bus, and stores the data into the 
respective Bus object. 

3.2.4 ASSIGNJTMETABLEO 

Once the required number of Stop objects, and Bus objects are 
created, the next step is to assign the time table to each Stop object, and Bus. object. 
This file opens the input file No. 3 and reads the time table for each Stop of the ring 
road network, and stores them in respective of Stops, and StopQue of Buses 
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by calling the getjimetab{ ) function of each Stop, and each Bus. The flow chart 
below explains the same process. 

3 . 2.5 OD_MANIPULATION( ) 

This function opens the input file No. 4, and reads the OD 
matrix of size MAX_STOP. The each element of OD matrix in the input file 
represents the daily demand from one Stop to the other. This function manipulates 
the OD matrix using normalization technique in such a way that each element 
represent the cumulative percentage of daily demand from one Stop to the other. The 
manipulated OD matrix will be used to assign the destinations randomly to the 
passengers by the function generate j)ass( ). 

3.2.6 DO__OPERATION( ) 

The functions explained so far in the previous sections open 
and read the four input files, create the Stop and Bus objects, assign the time table to 
them, and manipulate the OD matrix. But, they do no simulation, and scheduling. 
The function do_operation( ) is the important and most critical function, and it does 
the simulation and scheduling part of the modeling. This function contains sub 
functions, a few of which contain another level of functions. The sub functions used 
in this function are - 

1. find systemjime( ) 

2. generate _pass() 

3. get busJocations( ) 

4. findJbest_stop( ) 
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The function do_operation( ) will be explained with the help 
of flow charts in the following sections. Each of the sub functions will also be 
described. 

This function do_operation(^ ) first calls the function fi.xiA_systein ) to get 

system time ^ the time to start the simulation process. The system time, in this model 
is defined as the time 30 minutes before the minimum of first bus departure time from 
all Stops. Now let us try to study the functions used in od manipulationi , ). They 
are generate jpass{ ), fmd_best_stop ( ), and get_bus_locations { ). 

A) generate jpass(stop_no ) 

This function takes Stopjao, the id number of a Stop at which passengers are to 
be generated, and generates the passengers as per pre-decided probability 
distribution for a period of five minutes starting from system Jime, assigns the 
random destinations to he passengers, and stores them in PassQue of respective 

•'* I 

$ 

Stop.. The description of this function is given in section 2.5.4, which gives a 
precise over view of the function. However in this section let us study the same 
with the help of flow chart. The figure 3.3 shows the flow chart of the function 
generate j3ass( ). 

B) find best stopO 

The function ^m.d_best_stop { ) will be called at each increment of system time by 
five minutes, which evaluates the best strategy to provide an additional Bus on to 
the ring road. This function is explained in section 2.5.6 in detail. In this" section 
the flow chart is given to understand the same more clearly. The figure 3.4 shows 
the flow chart of the Unction findJ}est_stop { ). 
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C) get_bus_locations() 

This function at each increment of system time calculates the positions of Buses and 
plots the Buses on the ring road. The section 2.5.7 explains the function. The flow 
chart is presented in this section to enable the readers understand the logic used in 
this function 



\ / 


sort_stopJimes{ ) 
sort_bus_times ( ) 

( END 1 


Fig. 3.1 Flow Chart for Function assignjdmetable( ) 





START 
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fig. 3.2a Flow Chart for Function do_operation( ) ( . continued) 












51 







m 



Fig 3.3a Function to generate Gamma random variates 
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fig 3.4 Flow Chart for Function find_best_stop( ) 
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Fig. 3,5 Flow Chart for Function get_bus_location( ) 
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START 


set 1 = 1 


Stop[\\ store _prev_xy{ ) 
set 1 = 1 


1 += 1 
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|Yes 
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^ =0 . 


wtemjifne = pj 
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Bus X = next x 
Bus V = next v 


9 
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Bus x = next x 
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^0 
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time_diff = next_imie - p)evjime 
time Jag = system time - pi ev time 


next X = 


totjiist = prev y - next y 
move = tot_dist x time lag I time_diff 


tot_ciist = pi ev y - next y 
move = tot_ dist x time Jag / timejdrff 


next y = 
sPrevy y 

lYes 


Bus x = prevx"> 

Bus V = nrev v - move 

T ' 

I dra\\>_biis{) f END 


draw_bus( ) 


END 


lot_dist = V (prev x-next x)^ + (prev.y-next y)^) 
move = tot_dist x timejag / time_diff 
slope = ( prev y - next y ) / ( nprev.x - next y ) 

0 = tan''(slope) 



prevx< 

No \.next x/ 


Bus X = piev X - move x cos9 
Bus.y = prev.y - move x sin0 


Bus X = prev x + move x cos0 
Bus.y = piev y + move x sin0 


dra\vj)iis( ) 


END 


Fig. 3.6 Flow Chart for Function Jind bus _poisition( ) 
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CHAPTER IV 

MODEL APPLICATION 

- A CASE STUDY 


4.1 INTRODUCTION 


The model discussed in the pievious sections is tested foi a case 
study of New Delhi metiopolitan aiea The iing load netwoik of New Delhi city is 
consideied and the transit system on the iing road is simulated using this model The 
input for this model is extiacted fiom die data collected fiom piimaiy surveys, and 
secondaiy sources for the study ‘Optimal Design of Uiban Bus system’ conducted by 
RITES for New Delhi metiopolitan area (RITES, 1993). The data is available for the 
complete tianspoitation network of New Delhi metropolitan area in terms of link 
distances, daily demand fiom one node to the other for all 1200 nodes. The input thus 
extracted is fed to the model and the output is obtained. The following sections 
explain the extraction of input and generation of output both in graphics part and non- 


graphics part. 
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4.2 INPUT TO THE MODEL 

As explained befoie, the input to the model is exti acted fiom 
the data collected foi the study ‘ Optimal Design of Uiban Bus System ’ conducted 
foi New Delhi nietiopolitan aiea (RITES, 1993), This section explains how the 
lequned input is deiived fiom the complete data Fust all the bus stops on the iiiig 
load netwoik aie identified, the daily demand matrix foi those bus stops on the iing 
load is extiacted fiom that of the complete transportation netwoik, the major stops aie 
identified as terminals on the img road; the tiavel times on all the links of iing road 
network aie obtained, and a tentative lime table for the bus system on the iing load is 
generated 

4.2.1 IDENTIFICATION OF STOPS ON THE RINGROAD : 

Fust the img load is identified from the map of the complete 
netwoik of Delhi aiea, and then all the Stops on the ring toad ate noted. Foi Delhi 
netwoik with 1200 nodes the total numbei of nodes on the img road are found to be 
76. These 76 nodes aie veiy closely spaced and it would be complicated to display the 
movement of buses giaphicaliy on the scieen with such closely spaced nodes Flence 
some closely spaced nodes ate meiged ni to a single node, thus obtaining 31 nodes on 
the ring load network. These 31 nodes aie identified based on the link distances, and 
the daily demand from one node to the other. The model is then used to simulate the 
public tiansit system on the ting load netwoik of Delhi city with 31 bus stops on it 
The X and y coordinates of all the 31 stops aie obtained. The Table 4 1 gives the 
names of all 3 1 stops on the ring road network 
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Table 4.1 The major Stops on the Ring Road NetAvork of Delhi 


Stop Sequence No. 

Name of the Stop 

1 

RAJA GARDEN 

2 

PANJABI BAGH 

3 

ROTH AK ROAD 

4 

SHAKURPUR 

5 

WAZIRPUR 

6 

PITAMPURA 

7 

SHALIMAR BAGH 

8 

AZADPUR 

9 

MODEL TOWN 

10 

GURUTEJ BAGH NAGAR 

11 

KHALSA COLLEGE 

12 

KHYBERPASS 

13 

OLD SECRETERIATE 

14 

ISBT 

15 

YAMLINA BAZAR 

16 

RED FORT 

17 

1 1 P STADIUM 

18 

PRAGATI MAIDAN 

19 

NIZAMUDDIN 

20 

1 KALEKHA 

21 

BHAGA WAN NAGAR 

22 

ASHRAM 

23 

LAJP AT NAGAR 

24 

MOOL CHAND 

25 

SOUTH EXTENSION 

26 

A LI M S 

27 

BHII<LAJI 

28 

MOTI BAGH 

29 

DHAULA KUAN 

30 

BRAR SQUARE 

31 

MAYA PURI 


The Table 4.2 gives the information of stops with lespect to the 
desired piobability distribution of passengers aiiival at stops and its parameteis In 
this Table EX lefeis to Exponential distiibution , NR refeis to Noimal Distiibution, 
GM lefers to Gamma distiibution, In Table 4 2, par_l, and par_2 aie respectively fust 
and second parameters of the distribution. For Exponential distribution only pai_l is 

considered as p, and there will be no par_2. For Normal distribution par_l is 

« *2 

r» o ii nniacivt ftnrJ rkar 0 1 G pnncMprpH pf fhp deviatioTi. Fni 
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Gamma distiibulion par_l is consideied as a, the shape paiametei, and pai_2 is 
consideied as (3, the scale paiameter. 

Table 4.2 The Dcsiicd Probability Distributions of Passenger Arrivals at 
Stops 


Stop Sequence 
No, 

Desired 

Distribution 

par i of the 
Distribution 

par_2 of the 
Distribution 

1 

EX 

2.5 


2 

EX 

32 


3 

NR 

2 

0 5 

4 

EX 

3 


5 

GM 

1.2 

2 

6 

EX 

3.6 


7 

NR 

3 

05 

8 

EX 

2.7 



EX 

28 



EX 

3 1 


11 

NR 

25 

0.4 

12 

EX 

3 


13 

EX 

3 5 


14 

EX 

29 


15 

NR 

4 

1.33 

16 

EX 

39 


17 

EX 

27 


18 

GM 

1.0 

1 6 


EX 

3 1 



EX 

2.3 


21 

NR 

24 

.6 

22 

EX 

2.9 


23 

EX 

2.87 


24 

GM 

2 1 

1.3 

25 

EX 

2.1 


26 

EX 

28 


27 

NR 

3 

8 

28 

EX 

2 


29 

EX 

2 95 



EX 

3 8 


31 

NR 

2.7 

.5 
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4.2.2 FORMATION OF OD~MATRIX : Fiom the 

complete daily demand matux of size 1200x1200 of Delhi aiea first the daily demand 
matiix foi all 71 nodes on the iing load network is extracted. The OD matiix for the 
3 1 stops is then obtained by meiging the OD of smallei nodes (the nodes with less 
demand and closely spaced) with that of pie-identified stops Thus obtained OD 
matiix of size 31x31 foi all stops on the iing road netwoik of Delhi aiea Now each 


element of this matrix lepiesents the daily demand fiom one stop to the other, only on 
the ling load netwoik. The table 4 3 gives the total productions of all the Stops in 
teims of daily demand 

Table 4.3 Total Productions of all Stops 


Stop Sequence 
No. 

Total Productions 
of the Stop 

Stop Sequence 

No. 

Total Productions 
of the Stop 

1 

3560 

16 

716 

2 

4331 

17 

1593 

3 

1672 

18 

1645 

4 

2791 

19 

2141 

5 

1541 

20 

1129 

6 

1441 

21 

894 

7 


22 

487 

8 

3072 

23 

1045 

9 

2974 

24 

932 

10 

2613 

25 

1464 

11 

3573 

26 

837 

12 

3727 

27 

2260 

13 

3028 

28 

2781 

14 

383 

29 

1383 

15 

1339 

30 

1324 



31 

1 

1950 


4.2.3 IDENTIFICATION OF TERMINALS : 

After the stops on the ring road are identified, it is necessaiy to 
have few teiminals which could supply buses to the ting load. A terminal m this 
model is defined a stop which supplies the buses to the route network under 
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consiclcialion Each of the teiiiiinal will have some niimbei of available buses, which 
IS lefeiied to as tot_bus in the model The stops which aie not teiminals will not have 
any buses , i e, tot_bus foi such stop will be zero. With the initial time table given to 
the model as input, some numbei of buses will be scheduled from the each teiminal 
Consicleimg the lemaining buses at each teiminal, the model finds the best teimmal 
fiom which the additional bus could be scheduled 

Six leiminals aie identified out of 31 stops on the ring load 
based on the total demand fiom each stop to all the otliei stops It is assumed in this 
model that these teiminals supply buses to the ling load nelwoik, i e all the buses 
oiigiiiatc fiom any of the teiminals, move on the iing load seiving all the inteimediate 
stops and they finally stop at the lespective teiminals accoidmg to the pi e-detei mined 
timetable. The Table 4 4 gives the names of the selected teiminals on the iing road, 
and the number of buses available at theses teiminals 

Table 4,4 The Terminals on the Ring Road Network of Delhi. 


Terminal 

No. 

Name of tJic 

Terminal 

No. of Buses 
at Terminal 

1 

RAJA GARDEN 

18 

8 

AZADPUR 

8 

14 

I.S.B.T 

8 

18 

PRAGATI MAID AN 

8 

23 

LAJPATNAGAR 

8 

29 

DHAULAKUAN 

8 


4.2.4 ms INFORMATION : 

The simulation is initially started with 24 buses, foui buses 
from each of the six teiminals. Thus the FLEET__SIZE is equal to 24 before stalling 
the simulation The maximum value of the FLEET_SIZE is MAX_BUS, the total 
number of buses in the system, which is equal to the 43 in this case study. This inodel 
suggests a few more additional buses which will be added to tlie FLEET_SIZE, The 
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headway of these buses is kept as 20 minutes at each teiminal. The movement of 
buses on both diiections (Clock Wise, and Anti Clock Wise) on the ring load is 
consideied The Table 4 5 gives the infoiination of buses initially under opeiation is 
given in tei ms of its stalling teiminal, starting time, and its direction of movement 


Table 4.5 Information of Buses 
(CW - Clock Wise, ACW - Anh Clock Wise ) 


Bus 

Sequence No. 

Starting 

Terminal 

Starting 

Time (Hrs.) 

Direction of 
Movement 

1 

1 

610 

CW 

2 

1 

630 

ACW 

3 

1 

650 

CW 

4 

1 

710 

ACW 

5 

8 

605 

CW 

6 

8 

625 

ACW 

7 

8 

645 

CW 

8 

8 

705 

ACW 

9 

14 

605 

CW 

10 

14 

625 

ACW 

11 

14 

645 

CW 

12 

' 14 

705 

ACW 

13 

18 

600 

CW 

14 

18 

620 

ACW 

15 

18 

640 

CW 

16 

18 

700 

ACW 

17 

23 

615 

CW 

18 

23 

635 

ACW 

19 

23 

655 

CW 

20 

23 

715 

ACW 

21 

29 

600 

CW 

22 

29 

620 j 

ACW 

23 

29 

640 

CW 

24 

29 

700 1 

ACW 


4.2.5 PREPARING INITIAL TIMETABLE : 

An initial timetable is prepared to test the model, First tlie 
travel distance on each link is obtained from the data. For this case, iing load with 31 
stops on it, the total number of links aie 62, considering the movement of buses in 
both directions i.e clockwise, and anti clockwise. Assuming a constant speed of 15 
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kiiiph the tiavcl tune on each link are calculated jn teims of minutes. Having given in 
the input the staiting teiminal and diiection of movement of each of the 24 buses as 
shown in Table 4 5, the tentative timetable is piepaied The timetable piepared to test 
this model covcis the moining peak peiiod fiom 0700 to 1000 his Consideiing the 
diiection of movement of bus on img road, and its starting terminal the next 
destination of that bus is obtained, and the departure tune of the bus is computed 
using the following equation' 

Foi a bus, 

clepaituie tune at stopfi] = depaiture tune at stop[i -1] -f- tiavel time( i -1 , i ) 

+ boaidmg/alighting tune of passengeis at stop[i]. 

To develop initial timetable it is necessaiy to have the 
information of each link on the iing road netwoik, in teims of tiavel distance, and 
travel time Table 4 6 gives the details of each link on the iing road of Delhi city 
Tabic 4.6 Travel Distances, and Travel Times on Links of Ring Road. 


From - 
Stop No. 

To- 

Stop No. 

Link Distance 
(meters) 

Ti avcl time 
(minutes) 

Link Direction 

CW - Clock Wise 

ACW - AnLi Clock Wihc 

1 

2 

2400 

9 

CW 

2 

1 

2400 

9 

ACW 

2 

3 

700 

2 

CW 

3 

2 

700 

2 

ACW 

3 

4 

1270 

5 

CW 

4 

3 

1270 

5 1 

ACW 

4 

5 

1380 

5 

CW 

5 

4 

1380 ! 

5 1 

ACW 

5 

6 

780 

3 

CW 

6 

5 

780 i 

3 

ACW 

6 

7 

1050 

4 

1 

CW 

. 
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Tabic 4.6 cont’d 


From - 
Stop No. 

To- 

Stop No. 

Link Distance 
(meters) 

Travel time 
(minutes) 

Link Direction 

CW' Clock Wise 

ACW - Anil Clock Wise 

7 

6 

1050 

4 

ACW 

7 

8 

860 

3 

CW 

8 

7 

860 

3 

ACW 

8 

9 

2230 

8 

CW 

9 

8 

2230 

8 

ACW 

9 

10 

1380 

5 

CW 

10 

9 

1380 

5 

ACW 

10 

11 

680 

2 

CW 

n 

10 

680 

2 

ACW 

11 

12 

1100 

4 

CW 

12 

11 

1100 

4 

ACW 

12 

13 

1150 

4 

CW 

13 

12 

1150 

4 

ACW 

13 

14 

1760 

7 

CW 

14 

13 

1760 

7 

ACW 

14 

15 

1230 

4 

CW 

15 

14 

1230 

4 

ACW 

15 

16 

880 

3 

CW 

16 

15 

880 

3 

ACW 

16 

17 

2000 

8 

CW 

17 

16 

2000 

8 

ACW 

17 

18 

1680 

6 j 

CW 

18 

17 

1680 

6 

ACW 

18 

19 

2880 

11 

CW 

19 

18 1 

2880 

11 

ACW 

19 

20 

' 1600 

6 

CW 

20 

19 

1600 

6 

ACW 

20 

21 

600 

2 

CW 

21 

20 

600 

2 

ACW 

21 

22 

1250 

5 

CW 

22 

21 

1250 

5 

ACW 

22 

23 

1600 

6 

CW 

23 

22 

1600 

6 

ACW 

23 

24 

1200 

4 

CW 

24 

23 

1200 

4 

ACW 

24 

25 

1620 

6 

CW 

25 

24 

1620 

6 

ACW 

25 

26 

970 

3 

CW 

26 

25 

970 

3 

ACW 

26 

27 

1870 

7 

CW 

27 

26 

1870 

7 

ACW 

27 

28 

1710 

6 

CW 

28 

27 

1 1710 

6 

1 ACW 
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Table 4,6 cont'd. 


From - 
Slop No. 

To - Stop 
No. 

Link Distance 
(meters) 

Travel time 
(minutes) 

Link Direction 

CW- Clock Wise 

ACW - Anh Clock Wise 

28 

29 

2250 

9 

CW 

29 

28 

2250 

9 

ACW 

29 

30 

1590 

6 

CW 

30 

29 

1590 

6 

ACW 

30 

31 

3450 

13 

CW 

31 

30 

3450 

13 

ACW 

31 

1 

1510 

6 

CW 

1 

31 

1510 

6 

ACW 


4,2. 6 STORING THE RELEVANT DATA IN INPUT FILES : 


After pieparing the input data it is lequiied to store them Jii 
appropriate input files As mentioned in section 2 4 the input data must be stoied in 
five input files Foi this case study the input files contain the following data - 

1, Input file No. 1, ‘STOP_XY.IN’ contains the x and y cooidinates, the desiied 
cUstiibution and its paiameteis, and the total number of buses available foi each of 
the 31 stops on the iing road. The fust element of this file is 31, the number of 
stops undei consideiation The maximum numbei of buses available rn the system 
is equal to the sum of total number of buses of all stops This is taken as 48 in this 
case study (Refer Table 4 4) 

2. The input file No. 2, ‘BUS_lNFO.IN’ contains tlie starting terminal and diicction 


of movement of each bus. The fust element of this file is 50, which is the capacity 
of the bus taken in this model. In this file the infoimation of 24 buses is stoied, i e. 
the numbei of buses being opeiated on the ling load network before starting the 
simulation is 24. Thus the fleet size initially is equal to 24 The maximum value of 
the fleet size is the maximum number of buses in the system, i.e 48. 

3, The input file No. 3, ‘TIME_TAB.IN’ contains the initial timetable of the public 
transit system developed for this model. The format of this file is first stop 



numbei, then the numbei of depaitmg bus, and its depaituie time The timetable 
genciated as discussed in pieviovis section is stoied in this file Foi this case study, 
this file contains initial timetable foi all the 24 buses seiving the 31 stops on the 
ling load netwoik 

4 Ihe input file No 4, ‘LINKTIME IN’ contains the tiavel time on each link of the 
ling load netwoik. Fust the oiigin stop, then the destination stop, the tiavel time in 
minutes, and then the diiection of that link is the foimat that is used in this file 
hor this case study the tiavel times of all 62 links of the iing road aie stoied in this 
file 

5 The input file No. 5, ‘OD_MAT.IN’ contains the daily demand matiix foi all stops 
on the ling load netwoik. This matiix for this case study is 31x31 symmetiic 
mail IX. 

4.3 SIMULATION RUNS 


The simulation piocess can be started by executing the 
piogiam, i.e. by typing ‘sim3000’ at DOS piompt. The mam window will be 
displayed whieh will have four push buttons. Esc button - to stop the simulation and 
to exit to DOS; Enter button - to continue the simulation; Ctrl+S button - to display 
the infoimation of all stops at any instant of time, and Ctrl+B button - to display the 
information of all buses at any instant of time Tlie graphic display of this model is 
user interactive, and hence user can piess any of the four buttons of his interest There 
will be a Time window which displays the time of simulation. The idea of how the 
model works is aheady discussed in the previous chapters. 

Figuie 4 I shows display of Main window of this simulation 

model. The four push buttons Esc key. Enter key, Ctrl+S key, and Ctrl+B key can be 

a nj=*rfr^rm<3 Hiffp.ip.nt an.tinns. THp. Tiiyip 
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window ill the mam scieen displays the time 0535 Hrs. just befoie starting the 
simulation The Ring Road netwoik of Delhi with 31 major stops on it can be 
obscived All the vehicles, lepiesented by gieen coloi icons aie at then lespective 
teimmals as the simulation has just started, and their depaituie time has not yet 
leached. 

Figuic 4 2 shows the window which will be displayed when the 
model suggests an additional bus to be scheduled The two buttons to accept the 
additional bus, and to icject the additional bus can be activated by Space bur key, and 
Esc key icspectivcly. 

Figtiie 4 3 shows the Stop window, which will be displayed 
when the C7/'/+iS' key in the Main window is activated This window piovides tlte 
sufficient inlbimation of all stops such as the numbet of passengeis waiting at a stop, 
the average and total waiting time of passengeis, the sequence number and depaitiire 
time of next bus, ... etc. 

Figuie 4.4 shows the Bus window, which will be displayed 
when the Clrh-B key in the Main window is activated. This window displays the 
infoimalion of all buses such as the number of passengeis occupied the bus, the 
direction of inovement of bus on the iing toad, the sequence numbei of, and depaituie 
lime at previous stop, and next stop 
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4.4 RESULTS AND DISCUSSION 


1 he simulation foi the case study of Delhi ling load netwoik is 
caiiiccl ovei as discussed in section 4 3 The simulation piocess is started at 0600 his 
and teiminalecl at 1 130 his. thus coveiing tJie moiuing {jeaic period 0700 his to WOO 
lus, and the residts obtained will be discussed in the tbiiowing sections 

4.4.1 ADDITIONAL BUSES : 


As slated in section 4,2 4, theie are 24 buses in the system that 
aie being operated initially This model has simulated the public tiansit system on the 
ling load of Delhi city foi the moining peak pciiod and has suggested the 19 
additional buses to ply on the iing loacl The details of the additional buses scheduled 
on the ring load is given in the following table These additional buses oiiginate at 
thcii lespcctive stalling terminal, move along the ung load seiving all the Stops on the 
ring load, and finally teiminate at the then oiigm stop, i.e stalling teiminal The 
Table 4.6 gives the stalling teiminal, stalling time, and the diiectioii of movement of 
19 additional buse.s. 

Table 4.7 The Additional Buses Suggested by the Model 


New Bus No. 

Starting 
Terminal No. 

Starting 
Tune (Hrs.) 

Direction of 
movement 

25 

29 

610 

CW 

26 

18 

615 

cw 

27 

14 

620 

CW 

28 

8 

635 

ACW 

29 

1 

645 

ACW 

30 

23 

650 

CW 

31 

14 1 

705 , 

ACW 

32 

18 

720 

CW 

33 

29 1 

725 1 

CW 

34 

8 

730 

ACW 

35 

23 ! 

745 

ACW 

36 1 

1 

750 

ACW 

37 ! 

8 1 

810 

CW 
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Tabic 4.7 Cont’d 


New Bus No. 

Starting 
Terminal No. 

Starting 
Time Hrs. 

Direction of 
movement 

38 

29 

815 

ACW 

39 

18 

820 

ACW 

40 

14 

840 

CW 

41 

18 

845 

CW 

42 

23 

850 

ACW 

43 

1 

855 

CW 


4.4.2 SAVING IN WAIT TIME 


In this section the save in the total wailing time of each stop is piesented The model 
has suggested 19 new buses in addition to the 24 initial buses The additional buses 
aie suggested based on the saving in total wait time, evaluating the best strategy as 
explained In section 2 5 5, Thus the scheduling of these 19 additional buses on the 
ling load has proved that the total waiting time at each stop is reduced It has been 
found that the tola! wailing time at a stop when all tlie additional buses, suggested by 
the model are accepted is less than the total waiting time at that stop when no 
additional bus is accepted This is piesented in Table 4.7 


Table 4.8 Total Saving in Passenger Wait Times 


Stop No. 

Total Wait Time 

Total Wait Time 

Saving ill Total 

Rejecting New Buses 

Accepting New Buses 

Wait Time 


( Minutes ) 

( Minutes ) 

( Minutes ) 

1 

2752 

1932 

820 

2 

2491 

2281 

210 

3 

3043 

2530 

513 

4 

3287 

2382 

905 

5 1 

2566 

2369 

197 

6 : 

3573 

3109 

464 

7 1 

3642 1 

2963 

679 

8 ! 

2798 1 

2141 

657 

9 

2649 

2544 

105 
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Tabic 4.8 Cont’d 




Total Wait Time 
Accepting New Buses 
( Minutes ) 

Saving in Total 
Wait Time 
( Minutes ) 


2702 

2180 

522 


3902 

3377 

525 

1^1 

2894 

2397 

497 


3343 

2976 

367 

14 

3758 

3294 

464 

15 

2992 

2472 

520 

16 

3632 

3157 

475 

17 

3425 

3051 

374 

18 

4535 

3929 

606 

19 

2857 

2317 

540 

20 

3168 

2496 

672 

21 

2654 

2183 

471 

22 

2938 

2585 

353 

23 

2379 

1902 

477 

24 


2461 

596 

25 

2468 

1836 

632 

26 

3648 

3172 

476 

27 

3746 

3282 

464 

28 

4316 

3442 

' 874 

29 

5674 

4261 

1413 


3258 

2843 

,415 

31 

3494 

2657 

837 


4.4.3 NEW TIMETABLE 


The model has suggested 19 additional buses, and all these 19 
additional buses wcic accepted to ply on the ring road The saving in total wait time 
due to the scheduling of additional buses is discussed in the previous section. The 
model has prepared the new timetable which includes the tentative timetable with 24 
initial buses, and the timing of the 19 additional buses The buses with sequence 
numbcis from 1 to 24 are the buses in the initial fleet, and with sequence numbeis 
fiom 25 lo 43 are the newly scheduled buses. The revised timetable is presented in 


Appendix. 
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CHAPTER V 

CONCLUSIONS 


5.1 CONCLUSIONS 

This model has the capability for dynamic schedidmg of vehicles 
on a major coiiidoi of a meliopolitan city. This model simulates the public tiansit system 
on any route nctwoik, and suggests the optimal utilization of vehicles which minimizes 
the passengci waiting times at all stops. Consideiing the laige variations in passengei 
demand it is ncccssaty to have a transit system with dynamic scheduling of vehicles foi 
optimal design This model is flexible in terms of piobability distributions of passengei 
arrival at stops. Passengci can ariivc the stops as pei any one of Exponential, Noimal, and 
Gamma distributions, and the paranietei of these distiibutions aie as per user choice, 
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The model is developed with Object Oiiented Piogiammmg (OOP) 
appioach using C 1 he ntodel siippoits the usei iiiteiactive giaphics, so that the usei can. 
obtain the inroimation ol all stops and buses at any instant of time and/oi* quit the 
simulation at any instant of lime The locations of all buses on the route network can be 
displayed at every 5 minute interval of time by the concuirent mode of animation The 
softwaie is developed and luii on Tuibo C platfoim in an MS-DOS enviionment The 

size of the pioblem will be gieatly influenced by the computational capabilities, and the 
available memory 

This model is tested with a case study of Delhi city The img load 
network of Delhi mcliopolitaii aica is consideied, and the public tiansit on the ring load is 
simulated The data collected foi a study ‘Optimal Design of Uiban Bus System foi 
Delhi’ IS used paitly as the input foi this model A tentative timetable is prepaied, which 
is fed to the model. The model has levised the timetable, which includes the scheduling 
of additional buses suggested by the model. For the case study, the model has suggested 
19 new buses in addition to the 24 initial buses. The saving in wait time at all stops due to 
the schedule of these additional buses is also piesented 

5.2 LIMITATIONS 

This model has some limitations which can be overcome by minoi 
adjustments, and modilications to the softwaie The following limitations aie thought of 



• The niodcil can simulale the tjansjt operations on a majoi couidoi of city but not 
many coiiidois. 

• The model can not lake caie of ovci lapping loutes 

• The anlval latc of passcngeis at a slop is assumed to be constant in this model Wheie 
as most of the passengcis aiiive the stops based on the scheduled bus timings. 

• The model has not considered the stochastic iiatuie of bus auivals. 

• This model can display the locations of buses on the loute netwoilc at eveiy five 
minute inlcival, but not at any mleimediate time. 

5.3 SCOPE FOR FUTURE WORK 

Consideiing the limitations mentioned m the above section, theie is 

a lot of scope foi the fiiluic woik on this aiea 

• The model could be developed so as to nlake on-line decisions of optimal bus 
scheduling for a city loutc netAvoik. 

• The simulation for dynamic scheduling is Imiited to a single corridor in this model, 

which can be extended to many corridors with ovei lapping routes. 



APPENDIX 


Tinictablc revised by the model. 


Includes the schedules of additional buses. 


1 

2 

3 

1 

2 

3 

1 

2 

3 

Stop 

Bus 

Depaihiie 

Stop 

Dus 

Departure 

Stop 

Bus 

Departure 

No. 

No. 

Time ( IIis) 

No 

No. 

Time ( His) 

No. 

No. 

Time ( Hrs) 

1 

1 

610 

1 

5 

826 

2 

28 

657 

1 

21 

625 

1 

18 

827 

2 

3 

659 

1 

2 

630 

1 

11 

836 

2 

23 

714 

1 

25 

635 

1 

22 

847 

2 

10 

717 

1 

29 

645 

1 

32 

850 

2 

17 

724 

1 

3 

650 

1 

43 

855 

2 

8 

727 

1 

6 

656 

1 

1 

902 

2 

14 

733 

1 

23 

705 

1 

7 

906 

2 

13 

739 

1 

28 

706 

1 

20 

907 

2 

34 

752 

1 

4 

710 

1 

2 

922 

2 

26 

754 

1 

17 

715 

1 

24 

927 

2 

12 

757 

1 

10 

726 

1 

29 

937 

2 

31 

757 

1 

13 

730 

1 

35 

937 

2 

30 

759 

1 

8 

736 

1 

3 

942 

2 

33 

759 

1 

14 

1 742 

1 

39 

942 

2 

19 

804 

1 

26 

745 

1 

4 

1002 

2 

9 

805 


30 

750 

1 

41 

1015 

2 

16 

813 

1 

33 

750 

1 

37 

1031 

2 

18 

818 

1 

36 

750 

1 

40 

1031 

2 

15 

819 

1 

19 

755 

1 

36 

1042 

2 

27 

820 

1 

9 

756 

1 

38 

1042 

2 

5 

835 

1 

34 

801 

1 

42 

1042 

2 

22 

838 

1 

12 

806 

1 

43 

1147 

2 

11 

845 

1 

31 

806 

2 

1 

619 

2 

20 

858 

1 

15 

810 

2 

21 

634 

2 

32 

859 

1 

27 

811 

2 

25 

644 

2 

43 

904 

1 

16 

822 

2 

6 

647 

2 

2 

913 











1 


t 

2 

3 

1 

2 

11 

13 

816 

12 

20 

11 

20 

821 

12 

13 

11 

37 

825 

12 

37 

11 

26 

831 

12 

2 

11 

2 

836 

12 

26 

11 

30 

836 

12 

24 

11 

33 

836 

12 

30 

11 

19 

841 

12 

33 

11 

24 

841 

12 

19 

11 

9 

842 

12 

9 

11 

29 

851 

12 

29 

11 

35 

851 

12 

35 

11 

15 

856 

12 

39 

11 

39 

856 

12 

6 

11 

27 

857 

12 

15 

11 

6 

902 

12 

27 

11 

28 

912 

12 

28 

11 

4 

916 

12 

4 

11 

11 

922 

12 

11 

11 

32 

936 

12 

8 

11 

43 

941 

12 

32 

11 

8 

942 

12 

43 

11 

36 

956 

12 

36 

11 

38 

956 

12 

38 

11 

42 

956 

12 

42 

11 

34 

1007 

12 

34 

11 

41 

IlOI 

12 

41 

11 

40 

1117 

12 

40 

12 

5 

624 

13 

5 

12 

10 

636 

13 

10 

12 

14 

652 

13 

14 

12 

1 

700 

13 

1 

12 

7 

704 

13 

7 

12 

21 

715 

13 

12 

12 

12 

716 

13 

31 

12 

31 

716 

13 

21 

12 

25 

725 

13 

16 

12 

16 

732 

13 

25 

12 

18 1 

737 

13 

18 

12 

3 

740 

13 

3 

12 

23 ' 

755 

13 

22 

12 

22 

757 

13 

23 

12 

17 

805 

13 

17 


817 

13 

20 

813 

820 

13 

13 

824 

829 

13 

2 

828 

832 

13 

24 

833 

835 

13 

37 

833 

837 

13 

26 

839 

840 

13 

29 

843 

840 

13 

35 

843 

845 

13 

30 

844 

846 

13 

33 

844 

847 

13 

39 

848 

847 

13 

19 

849 

852 

13 

9 

850 

858 

13 

6 

854 

900 

13 

15 

904 

901 

13 

28 

904 

908 

13 

27 

905 

912 

13 

4 

908 

926 

13 

11 

930 

938 

13 

8 

934 

940 

13 

32 

944 

945 

13 

36 

948 

952 

13 

38 

948 

952 

13 

42 

948 

952 

13 

43 

949 

1003 

13 

34 

959 

1105 

13 

41 

1109 

1121 

13 

40 

1125 

628 

14 

9 

605 

632 

14 

27 

620 

648 

14 

10 

625 

704 

14 

5 

635 

708 

14 

14 

641 

712 

14 

11 ! 

645 

712 

14 

12 

705 

719 1 

14 

31 

705 

728 1 

14 

1 

711 

729 

14 

7 

715 

733 

14 

16 

721 

744 

14 

18 

726 

753 

14 

21 

726 

759 

14 

25 

736 

809 

14 

22 

746 




1 

2 

3 

1 

2 

14 

3 

751 

15 

18 

14 

20 

806 

15 

21 

14 

23 

806 

15 

25 

14 

17 

816 

15 

22 

14 

2 

821 

15 

3 

14 

24 

826 

15 

20 

14 

13 

831 

15 

23 

14 

29 

836 

15 

2 

14 

35 

836 

15 

17 

14 

37 

840 

15 

24 

14 

40 

840 

15 

29 

14 

39 

841 

15 

35 

14 

26 

846 

15 

13 

14 

6 

847 

15 

39 

14 

30 

851 

15 

6 

14 

33 

851 

15 

37 

14 

19 

856 

15 

40 

14 

9 

857 

15 

26 

14 

28 

857 

15 

28 

14 

4 

901 

15 

30 

14 

15 

911 

15 

33 

14 

27 

912 

15 

4 

14 

10 

917 

15 

19 

14 

8 

927 

15 

10 

14 

11 

937 

15 

15 

14 

36 

941 

15 

8 

14 

38 

941 

15 

36 

14 

42 

941 

15 

38 

14 

32 

951 

15 

42 

14 

34 

952 

15 

34 

14 

43 

956 

15 

12 

14 

12 

957 

15 

31 

14 

31 

957 

15 

32 

14 

41 

1116 

15 

43 

14 

40 

1132 

15 

41 

15 

9 

609 

16 

9 

15 

27 

624 

16 

27 

15 

14 

637 

16 

14 

15 

5 

639 

16 

5 

15 

11 

649 

16 

11 

15 

1 

715 

16 

16 

15 

16 

717 

16 

1 

15 

7 

719 

16 

18 


3 

1 

2 

3 

722 

16 

7 

722 

730 

16 

21 

733 

740 

16 

22 

739 

742 

16 

25 

743 

755 

16 

3 

758 

802 

16 

20 

759 

810 

16 

23 

813 

817 

16 

2 

814 

820 

16 

24 

819 

822 

16 

17 

823 

832 

16 

29 

829 

832 

16 

35 

829 

835 

16 

39 

834 

837 

16 

13 

838 

843 

16 

6 

840 

844 

16 

37 

847 

844 

16 

40 

847 

850 

16 

2‘8 

850 

853 

16 

26 

853 

855 

16 

4 

854 

855 

16 

30 

858 

857 

16 

33 

858 

900 

16 

19 

903 

913 

16 

10 

910 

915 

16 

15 

918 

923 

16 ■ 

8 

920 

937 

16 

36 

934 

937 

16 

38 

934 

937 

16 

42 

934 

948 

16 

34 

945 

953 

16 

12 

950 

953 

16 

31 

950 

955 

16 

32 

958 

1000 

16 

43 

1003 

1120 

16 

41 

1123 

612 

17 

9 , 

620 

627 

17 

14 

626 

634 

17 

27 

635 

642 

17 

5 

650 

652 

17 

11 

700 

714 

17 

16 

706 

718 

17 

18 

: 711 

719 

17 

1 

726 
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1 

2 

3 

1 

2 

3 

1 

2 

3 

25 

31 

856 

26 

3 

855 

27 ; 

25 

847 

25 

16 

912 

26 

16 

909 

27 


902 

25 

18 

917 

26 

23 

910 

27 

16 

902 

25 

41 

925 

26 

18 

914 

27 

18 

907 

25 

37 

941 

26 

41 

928 

27 

23 

917 

25 

40 

941 

26 

37 

944 

27 

41 

935 

25 

33 

952 

26 

40 

944 

27 

20 

947 

25 

20 

957 

26 

20 

954 

27 

37 

951 

25 

35 

1027 

26 

33 

955 

27 

40 

951 

25 

39 

1032 

26 

35 

1024 

27 

33 

1002 

25 

43 

1057 

26 

39 

1029 i 

27 

35 

1017 

25 

42 

1132 

26 

43 

1100 

27 

39 

1022 

26 

17 

628 

26 

42 

1129 

27 

43 

1107 

26 

22 

642 

27 

17 

635 

27 

42 

1122 

26 

13 

643 

27 

22 

635 

28 

22 

629 

26 

26 

658 

27 

13 

650 

28 

17 

641 

26 

30 

703 

27 

26 

705 

28 

13 

656 

26 

19 

708 

27 

2 

710 

28 

2 

704 

26 

9 

709 

27 

30 

710 

28 

24 

709 

26 

2 

717 

27 

19 

715 

28 

26 

711 

26 

24 

722 

27 

24 

715 

28 

30 

716 

26 

15 

723 

27 

9 

716 

28 

29 

719 

26 

27 

724 

27 

29 1 

725 

28 

19 

721 

26 

29 

732 

27 

15 

730 

28 

9 

722 

26 

5 

739 

27 

27 

731 

28 

6 

730 

26 

6 

743 

27 

6 

736 

28 

15 

736 

26 

U 

749 

27 

5 

746 

28 

27 

737 

26 

28 

753 

27 

28 

746 

28 

28 

740 1 

26 

4 

757 

27 

4 

750 

28 

4 

744 

26 

32 

803 

27 

11 

756 

28 

5 

752 

26 

10 

813 

27 

10 

806 

28 

10 

800 

26 

I 

815 

27 

32 

810 

28 

11 

802 

26 

7 

819 

27 

8 

816 

28 

8 

810 

26 

8 

823 

27 

1 

822 

28 

14 

816 

26 

14 

829 

27 

14 

822 

28 

32 

816 

26 

21 

830 


7 

826 

28 

36 

824 

26 

36 

837 

i 

27 

36 

830 

28 

38 

824 

26 

26 

26 

26 

26 

38 

25 

34 

12 

31 

837 

840 

848 

853 

853 

27 

27 

27 

27 

27 

38 

21 

34 

12 

31 

830 

837 

841 

846 

846 

28 

28 

28 

28 

28 

1 

7 

34 

12 

1 

828 

832 

835 

840 

1 840 




1 

2 

3 

1 

2 

3 

1 

2 

3 


28 

28 

28 

28 

28 

28 

28 

28 

28 

28 

28 

28 

28 

28 

28 

29 

29 

29 

29 

29 

29 

29 1 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

21 

25 

16 

18 

3 

23 

20 

41 

37 

40 

33 

35 

39 

43 

42 

21 

25 

22 

23 

17 

2 

24 

13 

29 

26 

6 

30 

33 

19 

9 

28 

4 

15 

27 

10 

5 

8 

14 

11 

36 

38 

32 

843 

853 

856 

901 

908 

923 

941 

941 

957 

957 

1008 

1011 

1016 

1113 

1116 

600 

610 

620 

640 

650 

655 

700 

705 

710 

720 

721 

725 

725 

730 

731 

731 

735 

745 

746 

751 

801 

801 

807 

811 

815 

815 

825 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

29 

30 

30 

30 

30 

30 

30 

30 

30 

1 30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

34 

12 

31 

1 

7 

16 

18 

21 

25 

22 

3 

20 

23 

41 

24 

35 

37 

40 

39 

33 

38 

42 

43 

21 

25 

23 

2 

17 

29 

13 

6 

28 

26 

4 

30 

33 

19 

9 

10 

15 

27 

8 

826 

831 

831 

837 

841 

847 

852 

852 

902 

912 

917 

932 

932 

950 

952 

1002 

1006 

1006 

1007 

1017 

1107 

1107 

1122 

606 

616 

646 

649 

656 

704 

711 

715 

725 

726 

729 

731 

731 

736 

737 

745 

751 

752 

755 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 

30 , 

30 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

31 

14 : 

5 1 

36 

U 

34 

12 

31 

32 

16 

1 

18 

7 

22 

3 

20 

24 

35 

41 

39 

37 

40 

38 

42 
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